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This paper gives a general treatment and proof of the direct conservation law method 
presented in Part I (see |^). In particular, the treatment here applies to finding the local 
conservation laws of any system of one or more partial differential equations expressed in 
a standard Cauchy-Kovalevskaya form. A summary of the general method and its effective 
computational implementation is also given. 



1 Introduction 

In this paper we present a general treatment of the direct conservation law method 
introduced in Part I (see Ref. |^). In particular, in Sec. || we show how to find the local 
conservation laws for any system of one or more PDEs expressed in a standard Cauchy- 
Kovalevskaya form. We specifically treat nth order scalar PDEs in Sec. ||. In Sec. |^ we 
summarize the general method and discuss its effective implementation in computational 
terms. 

In order to make the treatment uniform, it is convenient to work with Cauchy-Kovalev- 
skaya systems of PDEs as follows. 



Definition 1.1 A PDE system with any number of independent and dependent variables 
has Cauchy-Kovalevskaya form in terms of a given independent variable if the system is 
in solved form for a pure derivative of the dependent variables with respect to the given 
independent variable, and if all other derivatives of dependent variables in the system 
are of lower order with respect to that independent variable. 



Typically, scalar PDEs admit a Cauchy-Kovalevskaya form by singling out a derivative 
with respect to one independent variable, or by making a point transformation (more 
generally a contact transformation) on the independent variables. For example: the wave 
equation 

Utx = 

admits the Cauchy-Kovalevskaya form uu — u^x after the point transformation t — > t — x, 
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a; — > .T + i; the harmonic equation 

admits the Cauchy-Kovalevskaya form Uyy — —u^x with respect to y. A less trivial 
example is the Kadomtsev-Petviashvili equation [ p^ 

This equation admits two obvious Cauchy-Kovalevskaya forms: Uyy ~ T(wta: + {''^''^x)x + 
Uxxxx) which is a second-order PDE with respect to y\ and u^xxx — T^yy — utx — {uux)x 
which is a fourth-order PDE with respect to x. 

As examples which are more involved, consider the modified Benjamin-Bona-Mahoney 
equation Q 

Ut + {1+ u'^)ux - U.^xt = 0, 
and the symmetric regularized long wave equation [T^ ] 

Utt + Uxx + UUtx + UxUt + Uttxx = 0. 

As it stands the modified Benjamin-Bona-Mahoney equation is not of Cauchy-Kovalev- 
skaya form with respect to either t or x, since the t derivatives of u appear in both pure 
and mixed derivative terms, while the highest order x derivative of u appears in a mixed 
derivative involving t and hence is not in solved form. Nevertheless, if one makes the point 
transformation t ^ t, x x — t, then the modified Benjamin-Bona-Mahoney equation 
becomes Uxxx — Uxxt + u'^u.^ + ut = which now is of third-order Cauchy-Kovalevskaya 
form with respect to x. The situation for the symmetric regularized long wave equation 
is similar. It is not of Cauchy-Kovalevskaya form as it stands, but after one makes the 
point transformation t ^ t — x, x —^ x + t it is of fourth-order Cauchy-Kovalevskaya form 
with respect to t or x: Uutt + u^xxx - '2uuxx + (2 - u)utt -I- (2 -I- u)uxx + ut' - u^^ = 0. 

Many PDE systems can be handled similarly to scalar PDEs. For example, the vector 
nonlinear Schroedinger equation 

iut + Uxx ± f{\u\)u = 0, u = (u\ . . . , u") 

admits the first-order Cauchy-Kovalevskaya form ut = iuxx i if{\u\)u with respect to t, 
as well as the second-order Cauchy-Kovalevskaya form Uxx = —iut T f{\^\)u with respect 
to X. A less obvious example is Navier's equations of isotropic elasticity, 

KUxX + l^Uyy + {K- ^)Vxy = 0, 
{k - fJ,)Uxy + fJ-VxX + HVyy = 0, 

K — const, = const. This PDE system admits a second-order Cauchy-Kovalevskaya 
form with respect to x or y: Uxx = -f "ya + ~ '^)^xy and Vxx = -fvyy -f (1 - j^)uxy. 

In general any Cauchy-Kovalevskaya form of a system of one or more PDEs can be used 
with no loss of completeness in finding the conservation laws admitted by the system. 
Given a Cauchy-Kovalevskaya PDE system, we let t denote the independent variable in 
the derivative which appears in solved form in the PDEs, with the remaining independent 
variables denoted by "). In order to obtain the most effective formulation of 
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the direct conservation law method, it is convenient to express the system in its equivalent 
first-order (evolution) form with respect to t. 

Hence, we consider a first-order Cauchy-Kovalevskaya system of PDEs with N depen- 
dent variables u — (u^ , . . . ,u ) and n + 1 independent variables {t,x), 

du' 

= — +/(t,a;,«,9,«,...,a» = 0, a = l,...,iV (1.1) 

with X derivatives of u up to some order m. We use dxU, d^u, etc. to denote all derivatives 
of u"^ of a given order with respect to a:*. We denote partial derivatives d/dt and d/dx'^ 

by subscripts t and i respectively. Corresponding total derivatives are denoted by and 
D-. We let (^g)p denote the hnearization operator of g"^ defined by 

and we let {^g)^ denote the adjoint operator defined by 

oq on On 

(^XW^ -Q^W^~ D,{^ + + i-irD^,....^ (^-^ W^), (1.3) 

acting on arbitrary functions V'',W^ respectively. 

Throughout we use the summation convention for repeated lower-case indices; we use 
an explicit summation sign where needed for summing over non-indices. 



2 General treatment 



We start by considering the determining equations for symmetries and adjoint symme- 
tries. Suppose X is the infinitesimal generator of a symmetry leaving invariant PDE 
system (LI). We denote Xu^ — rf , which satisfies 



(2.1) 



for all solutions u{t,x) of Eq. (O). This linearization of Eq. (LI) is the determining 
equation for symmetries (point-type as well as first-order and higher-order type | pT| ) 
rf (t, X, u, du, . . . , d^u) of the PDE system (p^), where d^u denotes all jth order deriva- 
tives of u with respect to all independent variables t, x. The adjoint of Eq. (2J) is given 

by 

= -D,Lo^ + iCX^^, a = l,....N (2.2) 



which is the determining equation for adjoint symmetries uj^{t, x, u, du, . . . , d^u) of the 



PDE system (1.1). In general, solutions of the adjoint symmetry equation (2.2) are not 



solutions of the symmetry equation (2.1), and there is no interpretation of adjoint sym- 
metries in terms of an infinitesimal generator leaving anything invariant. 

In order to solve the determining equations for 77*^ and uj^, one works on the space of 
solutions of the PDE system. This means we use the PDEs to eliminate in terms of 

etc. In particular, without loss of generality, we are free to let rj'^ and uj^ have 
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no dependence on and its differential consequences. Let 



(2.3) 



which is the total derivative with respect to t on the solution space of PDE system 

(In particular, = 2?^ when acting on all solutions u{t,x).) Then the determining 

equations explicitly become 



0= vy + ic,)X 



drf f dvj 
~dt ^ 



dif 
du'' 



dg° 
'du' 



dg 
du 



dg'^ 



■B^n", a=l,...,7V (2.4) 



for rf {t^ X, u, dxU, . . . , d^u)^ and 

0= -i?,^, + (/:;)>^ 



dt \ du'' ^ du 
dg'' , dg'' 

^d7''''-''\d^''^ 



dg' 



(2.5) 



for ijJ„{t, X, It, dxU, . . . , d^u). The solutions of Eqs. ( |2.4| ) and ( |2.5| ) yield all symmetries 
and adjoint symmetries up to any given order p. 
We now consider conservation laws. 



Definition 2.1 A local conservation law of PDE system (IT) is a divergence expression 
L»t$*(t, X, u,du,..., d'^u) + D,^^\t, X, u,du,..., d'^u) = (2.6) 
for all solutions u{t,x) of Eq. (1.1); <&* and <&' are called the conserved densities. 



The conservation equation (|2.6|) holds as an identity if, for all solutions u{t, x) of 



Eq. (1.1), 



(2.7) 



for some expressions 9^{t, x, u, du, . . . ,d'^ ^u), ip"^^ {t, x, u, du, . . . ,d^ ^u) with -0'^ = 
—tp^^. Such conservation laws are trivial. Only the nontrivial conservation laws of the 



PDE system (IT) are of interest 



Definition 2.2 A local conservation law (2.6) is nontrivial iff the conserved densities do 



not satisfy Eq. (2.7) 



Any nontrivial conserved densities that agree to within trivial conserved densities are 
regarded as defining the same nontrivial conservation law. There is further freedom in 
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the form of conserved densities since we are clearly free to replace — —g" in $* and 
on the solution space of PDE system (pTl]). Thus, without loss of generality we can 
consider $* and to depend only on t, x, u, and x derivatives of u. We refer to this as 
the normal form of the conservation law, 

D^^\t, X, M, d^u, d!^u) + D^^\t, X, M, d^u, dl^u) = (2.8) 



for all solutions u{t, x) of PDE system (1.1). In normal form, the freedom corresponding 
to trivial conserved densities is given by 

$* + D,^e\ ^ - V^e' + D^V'^' (2.9) 

where — —-0-'' do not depend on Ut and differential consequences. 

All nontrivial local conservation laws (in normal form) can be shown to arise from 
multipliers on the PDEs (1.1) as follows. We move off the solution space of Eq. (1.1) and 
let u{t, x) be an arbitrary function of t, x. 

Definition 2.3 Multipliers for PDE system ( |l . l|) are a set of expressions 
{A^(i, x, u, du, . . . , d''u), . . . , A^(t, x, u, du, . . . , d'^u)} 

satisfying 

{u^ + g^A^ ^ + (2.10) 

for some expressions $*(i, x, m, du, . . . , d'^u) and ^\t, x, u, du, . . . , d'^u) for all functions 
u{t, x). 

Given a conservation law (2.8), consider -Dj$* + D^^\ Clearly this expression must 



be proportional to + g"^ and its differential consequences in order to satisfy Eq. (2. 
The terms arise only from 

where (C^t)^ = {d<S>'/du'') + {d^'/du1)D^ + • ■ • + (^^V^^n-iJ^n ' ' ' denotes the 
linearization operator of To organize these terms we use the identities 

= « + g'')K'{<^') - (£$0.3" + ^.r* (2.12) 

where F* is given by an expression proportional to u^+g'^ (and differential consequences), 
and where 

Eu" = - D^du^ + D^D^du^^ +■■■ (2.13) 
is a restricted Euler operator. Thus, we have 

= 9i<f * - (£$*)./ + ^.r* + « + (2.14) 

In order for the conservation equation ( ^ ) to hold, the terms 9t$* — {C^t)^g'^ which do 
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not involve + g'^ must cancel and therefore we have 

Then combining expressions ( PH ) and ( pl5| ) we obtain 



with 



(2.15) 

(2.16) 
(2.17) 



When u{t,x) is restricted to the solution space of PDE system (1.1), then vanishes 
and the divergence expression (2.16) reduces to the conservation equation (2. 



Hence, the expressions {£'„<t($ )} define multipliers {A^} yielding a conservation law 
( p.8[ ). Furthermore, since does not depend on Ut and its differential consequences, 
we see that each multiplier expression A^ is a function only of t,x,u, and x derivatives 
of u. Most important, these expressions A^ are invariant under a change in by a 
trivial conserved density ( ^.9[ ) since £'„<^ annihilates divergences I?j6'' where 9^ depends 

on t,x,u and x derivatives of u. (In particular, if in normal form is trivial, then A^ 
is identically zero, and conversely.) Thus we have the following result. 



Theorem 2.4 For the Cauchy-Kovalevskaya PDE system every nontrivial con- 

servation law in normal form ^.^ ) is uniquely characterized by a set of multipliers {A^} 
with no dependence on Ut and differential consequences, satisfying the relations ( 2.1(\) 
and ( 2.17 ) holding for all functions u{t,x). 

From this result it is natural to define the order of a conservation law (^^) as the order 
of the highest x derivative of u in its multipliers ( ^.17 ) . 

Theorem 2.4 is the starting point for an effective approach to find conservation laws 
of PDE system (1.1) by use of multipliers. The standard determining condition ||ll[] for 
multiplier expressions A^{t,x,u,dxU, . . . ,d^u) arises from the definition (2.10) by the 
well-known result that divergence expressions are characterized by annihilation under 
the full Euler operator 



This yields (by a straightforward calculation) 
= E^. «A^ + /A J = -D,A^ + {C;rA + {C^^M + /), a 



(2.18) 
.,N (2.19) 



where (£^)^p is the adjoint operator of the linearization operator {Cjy)^^ defined by 



and 



9A, 



dA^ 



du' 



dA„ 



du' 



A'ap 



D. ■■■D,. V 

tl 'tp 



.5A, 



dA 

(-l)PD. - D (— ^ 



(2.20) 



W) (2.21) 
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acting on arbitrary functions V'' , W'^ . Here the determining condition ( ^.19 ) is required 

to hold for all functions u{t,x), i.e., this is necessary and sufhcient for u^A^ + g'^A^ to 
be a divergence expression. We give a simple direct proof in Sec. [2.5| . 

We now show how to convert the determining condition for into a system of deter- 
mining equations that allow one to work entirely on the space of solutions of PDE system 
(1.1) to find A^. Furthermore, we show that the resulting determining system consists 
of the adjoint symmetry determining equation (2.5;) augmented by extra determining 
equations giving necessary and sufficient conditions for an adjoint symmetry to be a set 
of multipliers yielding a conservation law. 



2.1 Conservation law determining system 

In the determining condition (2.1£) for A^(t, a;, m, dxU, . . . , d^u) of order p consider the 
terms involving . These terms arise just from D^A^ and {C*^)^pU^, and so it follows 
that Eq. (2.19) is a linear polynomial in and differential consequences of with 



respect to x. Since is required to be an arbitrary function of t and x, Eq. (2.19) splits 
into separate equations given by the coefficients of , u^^, etc. It is convenient to organize 

this splitting in terms of + = G" and differential consequences + D^g'^ = Dfi'^, 
etc., which we refer to as the leading terms (all other terms in the splitting are then 
referred to as non-leading). Then the leading and non- leading terms in the splitting must 
vanish separately. 

To carry out the splitting of D^A^, we use the identity 

D, = v^ + « + /)a„P + « + Dy)d^. + ■■■ 



which yields D,A^ - V,A^ + {C^n". 



Consequently, the non- leading terms in Eq. ( 2.19 ) are given by 

0= -V,A^ + {CXA^ 
5A, 



dt 



(7-1,. 



g'cr p 

dA 



dw 
dg" 



9A„ 



D - D / 



/5/ 



This is the adjoint symmetry equation (2.5) with uj^ — A^ 



The leading terms in Eq. ( 2.19 ) are given by 

0=-(^A).pG'' + (/:;).pG"', a^l, 



dg' 



dU: 



■A. 



(2.22) 



(2.23) 



which we call the adjoint invariance condition on A^. Now since u"^ is required to be an 
arbitrary function of t and x, we observe that Eq. (2.23) splits into separate equations 
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given by the coefficients of G" , D-G'^ ,. . -,0-^ ■ ■ ■ D-^G"^: 
^1 dA dA„ 



wfiere C 



dA^ 



^^^„^,+i„ 9A^ 



+ (-l)P-«C^D, -D 



9A, 



, q^l,...,p-l 



11- --ip 



= - 



aA 9A„ aA„ 



V f, ]_ Vjy 

l,...,7V;p= 1,...,7V 

This establishes the following important splitting result. 



q q\(r-q)\ 



(2.24) 



Lemma 2.5 For A^ with no dependence on Ut and differential consequences, the Euler 
operator equation (2.1i) is equivalent to the split system of equations {2.2i} and ( 2.2j\ ), 
which are required to hold for all functions u{t, x). 



Consequently, by combining Lemma 2^ and Theorem 2A, we see that Eqs. ( 2.22 ) 
and ( 2.24 ) constitute a necessary and sufficient determining system for finding multipliers 
{A^}. The number of equations in this system is ^ J"^-^,^-*' + ) '■'^"j"^-'' . 



Theorem 2.6 For the Cauchy-Kovalevskaya PDE system (1-1), the multipliers for all 



nontrivial conservation laws in normal form (2.S) up to any given order p are the solu 



Hons A^(t, X, u, dxU, . . . , d^u) of the determining system consisting of the adjoint sym- 
metry determining equation {2.2i) augmented by the extra determining equations ( 2.2^ . 
In particular, Eq. (2.24) gives necessary and sufficient conditions for an adjoint symme- 



try to be a set of multipliers. 

In deriving the determining system for A^, we have eliminated Ut and its differential 
consequences. As a result, one is able to work equivalently on the space of solutions of 
the PDE system (^) in order to solve the determining system to find A^ . In particular, 
the same algorithmic procedures which one uses to solve determining equations for sym- 
metries can be used to solve the determining system for multipliers. Moreover, there is 
freedom in mixing the order of solving the determining equations in this system. A direct 
(naive) approach is to solve the adjoint symmetry determining equation first, then check 
which of these adjoint symmetries satisfy the extra determining equations. As illustrated 
in the examples in Part I, a more effective approach is to use the extra determining 
equations first. 

Remarks on the extra determining equations: 



There is a simple interpretation of the extra determining equations (2.24). From rela- 



tion (2.17) between multipliers and conserved densities, we observe that A^ is a varia- 



tional expression (i.e. it arises as an Euler-Lagrange expression from $ ). The well-known 
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necessary and sufficient (Helmholtz) conditions for an expression to be variational 
are that its linearization operator is self-adjoint, and thus is a variational expression 
if and only if it satisfies |0, || 

(^a).p = (^I).p, a,p^l,...,N. (2.25) 



The operator equation (2.25) is a linear polynomial in of degree p. We easily find that 



if it is decomposed into separate equations given by the coefficients of the polynomial, 



then the resulting equations are the same as the determining equations (2.24) 



Corollary 2.7 Multipliers for any first- order C'auchy-Kovalevskaya PDE system are 
completely characterized as adjoint symmetries with a variational form. 



Moreover, it is interesting to note that the determining equations (2.24) take the same 



form regardless of g"^ for all first-order Cauchy-Kovalevskaya PDE systems ( |l.l| ) 



2.2 Conservation lav^r construction formula 

We now give an integral formula that constructs the conserved densities <i>* and $' for 



any nontrivial conservation law in normal form ( |2.8|) in terms of its multipliers {A^}. 
The formula makes use of the identities |0] 

WA^afpy' - ViC^W^ = D^S'[V, W;g], (2.26) 

W''{C^),,V' - V'iCD^^W^ = D^S^'[V,W;A], (2.27) 

where 

7n—lm—£—l J^n^ 

£=0 fc=o 



(2.28) 



^=0 fc=0 

(2.29) 

which are trilinear expressions derived by manipulation of the linearization operators 



and adjoint operators. (Note, the terms in Eqs. (2.28) and (2.29) with I — Q or k ~ Q are 
understood to involve no derivatives of V and W , respectively.) 
To set up the formula, we first let 

u^^^ = Xu" + (1 - X)u'' (2.30) 

where {u"} are any functions of t, x. This defines a one-parameter A family of functions 
with M^-^j = u' and u^q-j — u' . Then we let 

^p["(A)] ^ ^p{t,x,u^^yd^Uf^^y . . . ,dlu^^^), (2.31) 
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K{t,x) = ({ 



-/["(A)])Ap[«(A)] 



A=0 



(2.33) 



Theorem 2.8 For the Cauchy-Kovalevskaya PDE system the conserved densities 

of any nontrivial conservation law in normal form are given in terms of the multipliers 
by 



d\{u'' -u)K^[u^^^]+t dXK{Xt,Xx), 



/ dXy'K{\t,Xx) + 



dx(^S"[u-u,A[Uf^^^];g[u^^^] 



S'[u - u,g[u^^^] - Xg[u] + (1 - A)^*; A[u^^,]] 



(A)J 



(2.34) 



(2.35) 



In applying the construction formula (2.34) and (2.35), we must fix a choice for the 
functions {u"^}- If the expressions and g'^ are nonsingular for u"^ ~ 0, then we can 
choose — and this simplifies the integrals. Moreover, if u"^ = u'^ = satisfies the 
PDE system (lA), then the K integrals vanish. 

In the case when the expressions and g"^ are singular at u'' = (for some p = 

1, . . . , N), we must choose vf ^ Q such that the expressions A^[{i] and g'^{u] are nonsin- 
gular. It is sufficient to fix a simple choice of u'' such that the integrals converge. Any 
change in the choice of vf changes the conserved densities only by a trivial conserved 
density (|2j|) . 

A simple proof of Theorem 2.8 is given in Sec. 2.3. 



2.3 Proofs of Main Equations 



Recall that, for first-order Cauchy-Kovalevskaya PDE systems (1.1), the proof of the 



determining system (2.24) for conservation law multipliers in Theorem 2.6 reduces, by 



Lemma 2.5, to the determining condition (2.19) involving the Euler operator. To conclude 
this section, we present a simple, direct proof of this determining condition ( 2.19| ) together 
with the construction formula (2.34) and ( ^.35| ) for corresponding conserved densities in 
Theorem |2j. The proof of Eg. (|2.19|) is based on an identity for linearization of the 
multiplier equation (2.10). 
We let 



V) = - 1)"" + 

be a one-parameter family of functions with ^^x) — ^ 
with 9u(-_^-)/9A = u for any functions u {t^x). 



(2.36) 

being an arbitrary function, and 



Proposition 2.9 For any given expressions A^[m] = A^(i, a;, tt, dxU, . . . , d^u), = 
X, u, dxU, . . . , d^u) and l>' [tt] — $*(t, x, u, dxU, . . . , d^u), the following identities 
hold by direct calculation: 

(«) ^(KA)t+5"['"(A)])A.['"(A)]) 

= « + (^g[«,,,])>'')^-['^(A)] + Ka)* +/["(A)])('Ca[„^^^]).p^'' 
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+D 



(2.37) 



where S denotes the trilinear expressions given by Eqs. and 

(«) ^ (^t**[«(A)] + A.*'[«(A)]) = + DMi[^^jy) (2.38) 

where (C^^t)^ and (/^.j,)^ denote the linearization operators of and respectively. 



Proof of the multiplier determining condition and conserved density con- 
struction formula: 

Suppose <i>* , are conserved densities of a conservation law in normal form 



From Theorem 2.4 the multipliers for the conservation law are given by = i?„»(<i> ) 



satisfying the multiplier equation (^.10| ) with = = . 

Since the multiplier equation ( ^.10[ ) holds for all functions u'{t,x), it must hold for 
the one-parameter family u^^-^ . We now take the derivative of the resulting left-side and 
right-side expressions of Eq. (2.1C) with respect to A. By Proposition 2.9, on the left- 
side we obtain Eq. (2.37), while on the right-side we directly obtain Eq. ( ^.3g ). These 
expressions (2.37) and ( 2.38| ) are equal for all functions w'^(t, x) and therefore hold iff the 
terms multiplying v" vanish and the total derivative terms involving v"^ are separately 
equal (by considering the terms u^, From the terms multiplying v"^ we have 

= -D.AJu^,^] + iC;^^j:A^[u^,^] + iCl^^J.Mx,,+9>^,,]). (2.39) 

This reduces when A = 1 to Eq. ( 2.19| ) and hence {A^} is a solution of the determining 
condition ( 2.19| ). 

Conversely, suppose {A^} is a solution of the determining condition (2.19). Then, 
by combining the two identities in Proposition 2.9, we see A^ satisfies the linearized 
multiplier equation 



^ (Ka)* + 5^["(A)])A.[«(.)]) -D,^ ¥[u^,^] +D^^ $'[«(,^] (2.40) 



d 



d ^, 



dX ' ^^1' ' dX 
with d<i>'[u^^^]/dX and a$'[M(^^]/aA defined by 



(2.41) 



(^i[.,.j)>" = (^^[^'A[V)];9KA)]]+^1^'"(A)t+9[«(A)];AK,,]])-i?,0^^ 



"(A) 



(A)t 



■'(A)J 



^A)J 



(2.42) 



for some expressions 0^ , tp 



-^^^ . We now undo the linearization to obtain the multi- 



plier equation (2.10) by integrating with respect to A as follows. We set v"^ ^ u' — u 
and so 



U 



(A) 



X(u — u ) + u 



2.43) 
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Then we use the fundamental theorem of calculus to obtain 



[u^ + g"" [u])K^ [u] = [u] + [u] + {u^ + g" [u])A^ [u] - Z?,** [u] - [ii] (2.44) 



where 



-S^lu — ii, u. 



-g[u,.±A[u,.A]), 



(2.45) 



(2.46) 



to within trivial conserved densities. Since Eq. ( ^.44 ) holds for all u{t,x), while u{t,x) 
is fixed, we must have 



D^¥[u]+D^¥[u] = (u^ +g^[u])AJu] = K{t,x) 



It is then simple to check that Eq. (2.47) is satisfied identically by 

¥[u]^t [ dXK{Xt,Xx),¥[u]=x' [ dXK{Xt,Xx) 



(2.47) 



(2.48) 



Thus, we find from Eq. ( 2.44 ) that {A^} satisfies the multiplier equation ( 2.1C ), with 
conserved densities given by Eqs. ( p.45D to ( |2.48|) Hence, by Theorem are mul- 



tipliers for a conservation law in normal form (2.< 



To obtain the construction formula ( 2.34 ) and ( 2.35| ) for the conserved densities, we 
move onto the solution space of Eq. (IT) and substitute u'^x)t ~ ^^9'^[''A + (1 ^ ^)''^t 
Eqs. ( |2l5| ) and ( pl^ ). The expressions $ [u] and $' [u] directly reduce to the formula 



for and $' 



□ 



3 Treatment of Nth order scalar PDEs 

Here we exhibit the conservation law determining system and construction formula for 
scalar PDEs of any order with one dependent variable u and n + 1 independent variables 
). We work directly with the scalar PDE expressed in an iVth order 
Cauchy-Kovalevskaya form 

G = -Q^+ 9{t. X. u,du,..., 9"u) = (3.1) 

where in this section d'^u now denotes all derivatives of u of order q, excluding t derivatives 
of u of order q > N and their differential consequences (i.e. the PDE is written so that 
the t derivatives of u of highest order appear in solved form). 

Clearly, without loss of generality, for conservation laws we are free to eliminate TVth 
order t derivatives of u (and differential consequences) in considering conserved densities. 

Definition 3.1 A local conservation law in normal form for a Cauchy-Kovalevskaya 
scalar PDE {^) is a divergence expression 

D^^\t, X, u,du,..., d'^u) + D,^^\t, X, u,du,..., d^u) = (3.2) 
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holding for all solutions u{t, x) of Eq. (p.l| 



A conservation law ( |3.2| ) is trivial if it holds as an identity (2/7) for some expressions 
0^{t,x,u,du, . . . ,d''~'^u), {t,x,u,du, . . . ,d'^^^u) with V'*-' = — for all solutions 
u{t,x) of PDE {"iA). Only nontrivial conservation laws (3^) are of interest. 



All nontrivial conservation laws (3^) of PDE (O) can be shown to arise from multi- 
pliers on the PDE, similarly to Theorem 2A. We move off the solution space of Eq. (3.1) 
and let u{t,x) be an arbitrary function of t,x. We use the notation d^u — d'^u/dt^ for 
pure t derivatives of u, and Ui = du/dx"^ , Uij = d^u/dx'^dx^ , etc. for pure x derivatives 
of u, and diui = d'^^^u/dt'^dx'^, d^Uij = d'^^^u/dt'^dx^dx^ , etc. for mixed t,x derivatives 
of u, with OqU — u and d^Ui — Ui. 



Theorem 3.2 For the Cauchy-Kovalevskaya scalar PDE (3.1), every nontrivial conser- 
vation law (S.i) is uniquely characterized by a multiplier A with no dependence on d^u 



and differential consequences. The multiplier satisfies the relations 

(5f M + g)A = 7^,$* + D^{^' - r) 

and 

A = ^($*) 

holding for all functions u{t, x), where 



E = 



- D, 



d{d, 



d{d, 



(3.3) 
(3.4) 

(3.5) 



is a restricted Euler operator, and is given by an expression proportional to u + g 
and its differential consequences. 



From Eq. (3.4) one can show that A is invariant under a change in $ by a trivial 



conserved density (2.7). (In particular, if $ is trivial, then A is identically zero, and 



conversely.) Consequently, it is natural to define the order of a conservation law (3.2) as 
the order of the highest derivatives of u in its multiplier (3.4). 

It is straightforward to derive both the determining system for multipliers A and the 
construction formula for conserved densities in terms of A by applying the results in 



Sec. |2.l| and p.2| to the scalar PDE (3.1) written as a first-order Cauchy-Kovalevskaya 
system (which we carry out later). 

In order to display the determining equations explicitly, we introduce the A^ -|- 1 ex- 
pressions 



no = A, 



k=l 



dg 



A-D 



dg 



dg 



A) ), <7=1,...,A^ 



A +. 



(3.6) 
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where T>^ is the total derivative operator with respect to t on the solution space of the 
PDE (3.1) as defined by eliminating d^u = —g and all differential consequences. (In 

particular, V^u = dfU, T>^u — dfu, etc., and V^u = —g.) Note that, if the order of Oq 
with respect to x derivatives of u is p, the order of D,q is at most p + mq. 



Theorem 3.3 For the Cauchy-Kovalevskaya scalar PDE i\3.\ ), the multipliers for all 
nontrivial conservation laws up to any given order p are the solutions 

A(t, X, u, du, . . . , d^u) of the determining system 



f^AT = 



and 



d{dlu) ~ d{d^u) 



k=l 



d{d^u^,...ij 



p 



k=q+l 



kl 



-(-If 



q\{k-q)\ 



■D: 



d{dl 



0, 



where p' ^ p + mk, j = 0, 1, . . . , TV - 1; fc = 0, 1, . . . , TV - 1. 



(3.7) 



(3.8) 



In this system, Eq. (3.7) is the determining equation for the adjoint symmetries A 



uj{t, X, u, du, . . . , d^u) of order p of the PDE (3J), explicitly 



(3.9) 



The extra determining equations (3.8) are the necessary and sufficient conditions for an 



adjoint symmetry to be a conservation law multiplier. Since Eqs. ( p.7[ ) and (3.8) do not 
involve d^u or any of its differential consequences, one is able to work equivalently on 



the solution space of the PDE (3.1) in order to find the solutions A. 

In order now to display explicitly the construction formula for the conserved densities 
in terms of the multiplier A, we first define the trilinear expression 

N-l 



JV — 1 y 

S'[V,W;F]^ ^ (dIv[ 
j=0 ^ 



d{dlu,) 

dF 



d{dli 



W) 



dF 



W) 



(3.10) 



d{dlu,j,) ' d{dluu,j,) 

depending on arbitrary functions V,W, F. Next we let 

""(A) = + (1 — X)u 

where u is any function of t,x. This defines a one-parameter A family of functions with 



(3.11) 
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Uj-j^-) = u and Uj-q-j = u. Then we define 



^J^^CA)] = ^9(*'^''"(A)'^"(A)' ■ ■ ■ '^''"(A))' 9 = 0,1,..., TV -1 (3.12) 
5["(A)] = 9it,x,U(^^ydu^^^,...,dPu^^^), (3.13) 
Kit,x) = id^u + g[u])no[u], (3.14) 



using Eq. (3.6) for flq in terms of A. 



Theorem 3.4 for the Cauchy-Kovalevskaya scalar PDE (3.1), the conserved densities 



of any nontrivial conservation law (S.i) are given in terms of the multiplier A by 



pl N-l .1 

4>*=/ dXy^idfu - dlu)nj[u.^.]+t dXK(\t,Xx), 
Jo j^Q Jo 



(3.15) 



^' = x'J d\yK{\t,\x) + J d\(^S'[u-u,no[u^^^];g[ 

+S'[u - u,5[>)] - AffM + (1 - X)d^u;no[u^^^]]). (3.16) 



In applying the construction formula ( 3.15| ) and ( p.l(^ ), we fix the function u so that the 

expressions A[u] and g[u] are nonsingular. In particular, if A[0] and g[0] are nonsingular 
then we can choose u — 0, which significantly simplifies the integrals. Moreover, if m = 



u = satisfies the PDE (3J), then immediately the K integrals vanish. A change in the 
choice of u alters the conserved densities only by a trivial conserved density (2.7). 
Conversion to a first order Cauchy-Kovalevskaya system: 



We now outline the proof of Theorems 3.3 and p. 4 using Theorems E.q and 2.S. To 



begin we write the scalar PDE ( p.l| ) in first-order (evolution) form (1.1) with respect to 
t as follows: 



1 2 
U = U,U 



dtu, . . . 
■,9 



N 



N-l„ 



N N 

~u ,g 



G = Otu — u = 0, . . . , G — OfU — u = 0, G = Otu 



N 



(3.17) 
(3.18) 
g = 0. (3.19) 



Through Eqs. ( 3.17 ) to ( 3.19| ) there is a one-to-one correspondence between nontrivial 
conservation laws (3.2) of the scalar PDE (^^) and nontrivial conservation laws in nor- 



mal form (2_^) of the equivalent first-order PDE system ( 3.19 ). The relation between a 
multiplier A of a scalar PDE conservation law and a set of multipliers {A^, . . . , Aj^} of the 
corresponding PDE system conservation law can be obtained by considering the adjoint 
symmetry equations of the scalar PDE (3.1) and the PDE system ( 3.19| ). Straightfor- 
wardly, from Eqs. (3.18) and (2.22), we have 







= -T^An-, a 



N-q+1 



,N-1 



(3.20) 
(3.21) 
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where £* ^ is the adjoint operator of the linearization operator Cq^g defined by 



dg 



dg 



dg 



(3.22) 



By solving Eq. ( 3.2l| ) for A2,...,A^ in terms of A-^ and comparing Eq. ( |3.20 ) with 
Eq. ( |3.7|) , we directly see 



A^ = A = r^o, A^ 



N~l- 



This establishes an explicit correspondence between A and {Aj^, , 
diately to Theorems 3.3 and 3.4 from Theorems |2.6| and E^. 



(3.23) 

Ajy} leading imme- 



Remarks on the determining system and construction formula: 



Theorems 3.3 and pA can also be established directly from Theorem 3.2 without use of 



the results in Sec. and 2.2. The main step in the proof of Theorem B.S is a polynomial 



splitting result analogous to Lemma as follows. 

The determining condition for a multiplier A of order p for the scalar PDE ( |3.l|) arises 
from the relation ( |3.3| ) by the result that an expression is a divergence if and only if it is 
annihilated by the full Euler operator 

d „ d „ d d „2 d 

du 



E„ 



dui 



d{dtu) 



^ dui 



J d{dtuj) 



d{dfu) 



+ ■ 



(3.24) 



This can be shown (by a straightforward calculation ||rT|) to yield 



= i?.((af u)A + .gA) = i-D.^A + CIA + C\ {d^'u + g), 



(3.25) 



which is required to hold for all functions u{t,x) (not just solutions of Eq. ( |3.l| )). The 



determining condition (3.25) is a polynomial in u, 9, 



8''' 



and differential 



consequences with respect to x. Furthermore, the terms in this polynomial have weights 



up to N, where we assign weight 1 to u (and x derivatives of u), 2 to 
(and X derivatives of d^~^^u), etc., and we add the weights of products (and powers) 



of d^u^df u, etc.. Now, since u is required to be an arbitrary function of t and x, 
the polynomial splits into separate determining equations given by the coefficients of the 
various weight terms involving d^u, d^~^^u, . . . , d^^~^u (and differential consequences 
with respect to x). It is convenient to organize the splitting by working in terms of 
u + 5 = G, 5f + V^g = Vfi, 5f + D^g = Dfi, d^+\, + Dp^g = Dpfi, 



etc.. The terms of weight yield the adjoint symmetry determining equation (3.7) and 
the terms of weight 1 up to iV yield the extra determining equations (3.8) on A. This 
derivation is illustrated in the second example of Part I. 

The construction formula for conserved densities $* and of a conservation law for 
PDE (3.1) is obtained by inverting the Euler operator equation ( 3.25| ) as follows. Since 
Eq. (3.25) holds for arbitrary functions u(t,x), it must hold with u replaced by the 
one-parameter family Uf^^-^ = Au + (1 — X)u. This yields 



o = (-i?,)-A[.(,^] + /:; 



i9t u,^)+g[u,J). 



(3.26) 



We multiply Eq. ( |3.26D by u ~ u and then rearrange the terms which involve total 
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derivative operators coming from C* and £^ . This leads to the formula 



N-l 



+S'[u -u,d^Uf_^) +ff["(A)];^o[W(A)]]) = ^ ((5fw(A) +5["(a)])A["(a()?)27) 



Next we integrate from A = to A = 1 and apply the fundamental theorem of calculus. 
Using the identity D^(t dXK{Xt, Xx)^ + [x' dXX'^KjXt, Aa;)] = K, and finally 
moving onto the solution space of the P DE (|3.lD , we obtain the conservation law (3.2) 
with and given by Eqs. (|1|) and { ^.16\) . 
Remarks on variational principles: 



Definition 3.5 A Cauchy-Kovalevskaya scalar PDE (3.1) is called variational if it arises 
from an action 



S = J (^L{t,x,u, du,...,d^u)^dtdx 
by variation with respect to m, 

G^E^{L)^d^u+g. 



(3.28) 
(3.29) 



The well-known necessary and sufficient condition for existence of an action (3.28) 
is that 



D 



(3.30) 



i.e. N must be even and g must have a self-adjoint linearization. This condition is equiv- 



alent to requiring that the determining equation for symmetries of the PDE (3^) is 
self-adjoint. 



In the case when PDE ( p.l| ) is variational, Theorem 3.3 combined with Noether's 
theorem ||, |ll[] shows that the extra determining equations ( |3.8| ) constitute necessary 
and sufficient conditions for a symmetry of the PDE ( 3T ) to leave invariant the action 
( |j.28 ) to within a boundary term. In particular, if Xu — ri{t,x,u,du, . . . ,d^u) is a 

symmetry of order p, then XS = J {D^O* + D-9'^)dtdx holds for some expressions 9^ and 
9^ iff A = 7/ satisfies Eq. (3.8) and hence 77 is a multiplier yielding a conservation law 



(p^ of PDE (3.1) 



4 Summary and concluding remarks 



For any Cauchy-Kovalevskaya system G of one or more PDEs, Theorems 2.6, 2.8 and 
Theorems ^3, yield an effective computational method to obtain all local conservation 
laws (up to any specified order). The method is summarized as follows: 

1 . Linearize G to form its linearized system £, which is the determining system for the 
symmetries of G. 

2. Form the adjoint system £* of £, which is the determining system for the adjoint 
symmetries of G. 
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3. Form the extra system h comprising the necessary and sufficient determining equa- 
tions for an adjoint symmetry to be a muhipher for a conservation law of G. 

4. Solve the augmented system £*Uh. This is the determining system for the multipliers 
that yield all nontrivial local conservation laws of G. 

5. Use the explicit construction formula to obtain the conserved densities arising for 
each solution of the system £* U h. 

The linearized system of G is sclf-adjoint {£ = £*) if and only if G is variational, in 
which case solutions of £* are solutions of £. Then the extra system h is equivalent to 
the condition for symmetries to leave invariant the action for G. In general, if G is not 
variational then solutions of £* are not solutions of £. 

The systems £, £* , h, and £* U h are all linear overdetermined systems which are solved 
working entirely on the space of solutions of G (i.e. a leading derivative of the dependent 
variables in G is eliminated). There exist algorithmic procedures |^ to seek solutions of 
£. These procedures can be readily adapted for seeking solutions of h, and £* (J h. In 
general, £* \Jh is more overdetermined than £ and hence is typically easier to solve. More 
significantly, one can choose appropriate mixings of the determining equations in £* and 
h to solve £* U h effectively. 

One can also use specific ansatze to seek particular solutions of £*Uh, such as restricting 
the form of highest derivatives of the dependent variables of G allowed in the solution. 
For example, familiar conservation laws such as energy invariably arise from the simple 
ansatz of seeking multipliers restricted to be linear in first derivatives. 

In general it is important to note that solutions of £* are not necessarily solutions of 
h and hence £* does not determine a conservation law multiplier. This typically occurs 
for scaling symmetries of systems G in the case £* = £ (i.e. self-adjoint), and for point- 
type adjoint symmetries (first-order and linear in derivatives of dependent variables) of 
systems G in the case £* = —£ (i.e. skew-adjoint). Examples are utt — u^x -\- = Q 
which has u + tut + xux as a solution oi £ = £* but not a solution oi h\ ut + Uxxx = 
which has solution of £* = —£ but not a solution of h. Rcf. [|| exhibits several 

ODE examples in which nontrivial adjoint symmetries are not multipliers. The need for 
the extra conditions h to determine multipliers has not been clearly recognized in the 
literature (e.g. [^). 

The chief aspect of our method compared to other existing treatments of PDE con- 
servation laws (e.g. [0, [ll], H, HJ]) is the explicit delineation of the linear determining 
system £* U h which incorporates (and identifies) the necessary and sufficient conditions 
for adjoint symmetries to be multipliers, without moving off the space of solutions of 
the given PDE(s) G. Consequently, one can calculate multipliers of conservation laws by 
effective algorithmic procedures. Moreover there is the added computational advantage 
of allowing the determining equations in the adjoint system £* and the extra system h 
to be mingled to optimally solve the determining system £* U h, as illustrated by the 
conservation law classification results for the PDE examples in Part I. 
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